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Abstract 

We report the observation, by small-angle-neutron-scattering (SANS), of magnetic flux lines - 
"vortices" in super-clean KFe2As2 single crystals. The results show clear Bragg spots from a well 
ordered vortex lattice, for the first time in a FeAs superconductor. These measurements can give 
important information about the pairing state in this material, because the spatial variation of 
magnetic field in the vortex lattice reflects this pairing. With field parallel to the fourfold c-axis, 
nearly isotropic hexagonal packing of vortices was observed without VL-symmetry transitions up to 
high fields, indicating rather small anisotropy of the superconducting properties around this axis. 
This rules out gap nodes parallel to the c-axis, and thus d-wave and also anisotropic s- wave pairing. 
The strong temperature-dependence of the scattered intensity down toT « T c further indicates 
either widely different full gaps on different Fermi surface sheets, or nodal lines perpendicular to 
the axis. 

PACS numbers: 74.25.Uv, 74.70.Xa, 74.20.Rp, 74.25.-q 



Understanding the mechanism and symmetry of electron pairing in iron-based super- 
conductors represents an important challenge in condensed matter physics. Soon after the 
discovery of high temperature superconductivity in iron pnictide materials, many theorists 
suggested that pairing in these materials arises from a magnetic interaction [IJ. It is pro- 
posed that in optimally-doped materials, there is a repulsive interaction between holes on 
Fermi surface (FS) sheets at the center of the Brillouin zone and the almost equal number 
of electrons on other sheets at the corners. The simplest form of pairing that could arise 
from this interaction has the opposite sign of the pairing order parameter for the two types 
of carrier. Such a pairing state may be described as extended s- wave or s± since the order 
parameter has the same symmetry as the crystal, and is nonzero on all parts of the FS, 
despite changing sign. In this respect it may be contrasted with rf-wave pairing in cuprate 
materials, in which the order parameter has a lower symmetry than the underlying crystal, 
changes sign around a single sheet of FS and therefore has nodes. Other pairing states have 
been proposed for the pnictides j2j [3] which maintain s± symmetry, but in which the order 
parameter also changes sign on at least one FS sheet. Our material, KFe 2 As 2 , belongs to the 
'122' family of iron pnictide superconductors; the optimally-doped BaFe 2 (Asi_ x P x ) 2 member 
shows clear indications of nodes in the order parameter [H [5] which may imply either nodal 
s± pairing or alternatively rf-wave gap symmetry [3j. On the other hand, in optimally-doped 
Bai_ x K r Fe 2 As 2 with T c > 35 K, the order parameter appears to be fully-gapped on all parts 
of the FS, although with large and small gaps on different FS sheets [3 [7]. Optimally-doped 
BaFe 2 _ x Co x As 2 is similarly fully-gapped [8j. Such observations raise the question whether 
all the iron-based superconductors, even within the 122 family, have the same gap symmetry 
and pairing mechanism. We note, however that nodal gaps may be observed in the 'clean' 
materials with long mean free paths (which allow dHvA observations P ID]), and full gaps 
may be observed in 'dirtier' materials [6-8J. This indicates that low-energy states due to 
electron scattering are not a major cause of 'nodal' behavior. 

The vortex state in superconductors reflects the electron pairing state. However, SANS 
experiments up till now on pnictide superconductors have shown disordered VLs fTTl - fT3] . It 
is likely that nano-scale compositional inhomogeneity is a major source of vortex pinning 
and disorder in doped crystals. In KFe 2 As 2 , with Ba fully replaced by K, we have an over- 
doped but clean 122 material. Here we show how SANS measurements on the vortex lattice 
(VL) in this material can make a contribution to the debate on the gap structure in pnictide 



superconductors . 

For the present study, single crystals of stoichiometric KFe 2 As 2 were grown by a flux 
method. Detailed methods are described elsewhere [14] . SQUID magnetization measure- 
ments reveal that T c is 3.6 K with a 10 - 90% width of 0.2 K, and resistivity measurements 
give a typical residual resistance ratio <~ 400 [H]. These facts indicate the system is very 
clean and indeed dHvA signals are observed fTO] . 

SANS experiments were performed on instruments Dll and D22 at the ILL France and 
SANS-I at the SINQ facility in PSI, Switzerland. A mosaic of co-aligned crystals with the 
c-axis and one of the tetragonal [100] axes in a horizontal plane was glued to aluminium 
sample-plates with hydrogen- free CYTOP® varnish and cooled to a minimum of 1.5 K using 
a 4 He cryostat at Dll and to 50 mK using dilution refrigerators at D22 and SANS-I. The 
total mass of the crystals for the ILL experiments was about 100 mg but 250 mg at PSI. 
Using a cryomagnet, fields from 0.1 to 0.9 T were applied parallel to the crystal c-axis and 
approximately parallel to the incident neutron beam. Neutrons with a wavelength of X n = 
10 A were used and a position-sensitive detector was set at 15 m (ILL) and 18 m (PSI) 
from the sample position to obtain a reasonable q resolution. Except when measuring the 
main beam intensity, the undiffracted beam was caught on a neutron-absorbent beamstop. 
The sample had sufficiently weak pinning that the VL could be formed by applying the field 
well below T c ; field-cooling through T c was not required. Indeed, the VL perfection was 
slightly improved by oscillating the value of the field by ~ 0.01 T about its final value at 
base temperature, before taking measurements. To observe the VL diffraction, the sample 
and cryomagnet were rocked together either vertically or horizontally by small angles about 
the B || beam direction to bring the various diffraction spots through the Bragg condition. 
Backgrounds were taken with no VL present, either by heating above T c or by removing 
the applied field. Note that the real-space VL nearest-neighbor positions may be visualized 
by rotating the diffraction patterns by 90° about the field axis and adding another spot at 
the central main beam position. The SANS data were displayed and rocking curves were 
analyzed using the GRASP analysis package fT5] . 

Figures la-d show VL diffraction patterns with clear Bragg spots, measured at T = 50 
mK and 1.5 K at selected magnetic fields applied parallel to the fourfold [0 1] crystal 
axis and approximately parallel to the neutron beam. The 12 spots arise from a mixture 
of two domain orientations of hexagonal VLs, which have nearest-neighbor vortices along 



the tetragonal [1 0] and [0 1 0] axes, and are degenerate [Hj. Although the angles of the 
spots change with increasing field (Fig. le) and the observed intensity becomes weaker, the 
VL symmetry does not change for all conditions under which a VL signal is observed (see 
Fig. If) [16J. In this respect, KFe2As2 is unlike almost all conventional and unconventional 
superconductors with a fourfold axis [T7H2T]. In conventional materials, anisotropy of the 
Fermi velocity leads to VL phase transitions as the field is increased [22], so we conclude 
that the FS anisotropy is weak in KFe 2 As 2 . First-principles calculations using Eilenberger 
theory show that gap anisotropy has similar effects [23j [24] , and that a rf-wave nodal gap 
gives rise to a square VL at a small fraction (<~ 0.15 at T = 0.5T C ) of B c2 [25J. These 
calculations are confirmed by experiments on d-wave materials [201I2T] . Since in KFe2As2, a 
square VL is absent up to 0.9 T - approximately half the upper critical field B c2 at 50 mK, 
we rule out a nodal order parameter with the gap varying around the basal plane. 

Figure 2a shows the £>-field-dependence of the VL 'Form Factor' (F) derived from the 
integrated intensity [14] of the 'on axis' VL diffraction spots at 50 mK and at 1.5 K. F is a 
measure of the spatial variation of the field inside the mixed state; it generally decreases at 
large field as the vortex cores begin to overlap. These data were fitted to a modified London 
model with core/non-local corrections [26] : 

F = r ^exp(-c^ 2 ) (1) 

Here, A is the London penetration depth, £ is the coherence length and c is an empirical 
core cut- off parameter. The good fit suggests a conventional field-dependence without 
multi-gap effects [27] and allows us to extrapolate to obtain the zero-field value of A. The 
best fit at 50 mK gave A = 203 nm and c = 0.52, using £ = 13.5 nm (from B c2 = 1.8 T), 
and at 1.5 K: A = 240 nm and c = 0.55, with £ = 15.9 nm (from B c2 = 1.3 T). Random 
errors in A were small compared with probable systematic (calibration) errors ~ 5%. 

The low temperature value of A may be combined with normal state properties of 
KFe2As2 [2] to confirm that this material has a strongly enhanced carrier mass and that our 
samples are 'clean', i.e. the electron mean free path, £ 3> £ , the coherence length. Further- 
more, in Figure 2b, F for both on- and off- axis spots is plotted, but against q ) rather than 
field because the VL distortion gives the off-axis spots a larger q than the on-axis spots. We 
see that all results lie on nearly the same line, confirming the essential basal-plane isotropy 
of the pairing. 



Measurements of the intensity of the flux lattice diffraction signal at low fields are equiva- 
lent to measurements of the magnetic penetration depth and give its value and temperature- 
dependence. Hence, the temperature-dependence of the penetration depth (or equivalently 
superfluid density) is a direct measure of the effects of thermal excitation of quasiparticles 
over the gap. Figure 3a shows the temperature dependence of the peak intensity at our 
lowest field of B = 0.1 T. Figures 3b-d show fits of these data to different gap functions, 
(which are described in full in ref. [14]). From Fig. 3b, it is clear that a single full gap 
cannot fit the low temperature behavior. Also in Fig 3b, we find that a simple nodal gap, 
ignoring nonlocal effects [28] deviates from the data; however, a single nodal gap including 
these effects (Fig. 3c) can reproduce the data very well, with only one additional fitting 
parameter. Alternatively, we can obtain a good fit if at least three full gaps of very dif- 
ferent magnitudes are included, as shown in Fig. 3d. Although this uses a rather large 
number of parameters for the fit, we can not rule out this possibility because KFe 2 As 2 has 
multiple FS sheets jTQl [29j [30] , which could support multi- full-gap states. In any case, the 
low-temperature behavior is dominated by the minimum values that the gap has at low 
temperatures; the strong temperature-dependence of the intensity of our diffraction signal, 
continuing down to at least T <~ 0.02T C , indicates that KFe 2 As 2 has a range of gaps extend- 
ing down to very low values. Probably it has a nodal gap, which is also suggested by recent 
penetration depth [30j and other measurements [321 [33] performed down to various fractions 
ofT c . 

We now consider the fitted values of the ratio A(T = 0)/fcsT c . The "largest" gap in a 
material may give a larger ratio than the weak-coupling isotropic BCS value of 1.76: either 
because of gap anisotropy or strong- coupling. However, it is not possible for the largest 
gap in a superconductor to give a smaller ratio. The small values of the maximum full gaps 
in Figs. 3b & d argue against gaps without nodes. Even the nodal gap models give values 
of the ratio close to the weak-coupling limit. We conclude that unlike the optimally-doped 
compositions, KFe 2 As 2 is a fairly weakly- coupled superconductor. 

Finally we discuss the gap structure, bearing in mind that the temperature-dependence 
of the signal argues for a nodal state, but the VL structure measurements rule out vertical 
line nodes in the gap. Note that T c in the (Ba/ K)Fe 2 As 2 system appears to decrease 
continuously and monotonically as the doping is increased above optimal [34 j • This strongly 
suggests that the symmetry of the order parameter does not change with doping, even 
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though the prominent electron sheets in optimally-doped materials are replaced by small 
hole pockets in KFe 2 As 2 P2J E9j [30]. The multi-full-gap scenario may be consistent with 
this if essentially isotropic gaps with widely different values develop on each of the FS sheets. 
Alternatively, we propose an axially isotropic nodal pairing state to account for our results. 
In this model, the gap nodes circulate around the approximately cylindrical sheets of the FS, 
with the order parameter changing sign as a function of k Z) but not in the k x -k y plane. For 
such a state, the pairing interaction should be operative in the z-direction, which requires 
some three-dimensionality in the electronic structure; dHvA measurements [10J and the 
moderate anisotropy of H c2 [35j shows that some of the FS sheets in KFe 2 As 2 indeed have 
a three-dimensional character. We further propose that as the doping in the (Ba/K)Fe 2 As 2 
system is increased beyond optimal, the in-plane pairing interaction becomes weaker, and 
the nodeless gap changes into a /^-dependent nodal gap. These two states have the same s± 
symmetry, and the node positions are "accidental" (not symmetry-determined) [36] so the 
system may go continuously from one to the other as doping or carrier scattering are varied. 
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FIG. 1: Small angle neutron scattering patterns, VL structure and B-T phase diagram for 
KFe2As2- a-d: Diffraction patterns from the VL under the following conditions of applied field and 
temperature: (a) 0.1 T, 50 mK; (b) 0.2 T, 1.5 K; (c) 0.4 T, 1.5 K; (d) 0.9 T, 50 mK, obtained on 
instruments Dll and D22 at ILL and SANS-I at PSI. The diffracted intensity of a spot could be 
maximized by tilting the sample and field together onto the Bragg angle [14J. Backgrounds ob- 
tained without VL present were subtracted. To improve visibility, the data were smoothed with a 
Gaussian of width comparable to the instrument resolution, and Poisson noise near the main beam 
was masked. The higher field patterns are mosaics of results taken at the maximum intensities of 
four sets of three spots, and at 0.1 T a sum over tilts is shown. White circles in the figures indicate 
the distance from the main beam at which regular hexagonal VL spots would appear. The offsets 
of the spots from the circles indicate that the VL is somewhat distorted from a pure hexagonal 
structure, (e): Field dependence of the apical angle, /3, of the VL structure at T = 50 mK. (for 
definition of (3 see Figure la), (f): The B C 2-T phase diagram (from magnetization data with H || 
[0 1]) indicating all conditions under which VL measurements were carried out. 
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FIG. 2: (a): Field dependence of the Form Factor (F) of the 'on axis' VL diffraction spots at T = 
50 mK and 1.5 K. The theoretical fit is described in the text, (b): F for on- and off-axis spots at 
50 mK and 1.5 K (plotted versus q- values, which are different for the two types of spot). 
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FIG. 3: (a): Temperature dependence of peak intensity at B = 0.1 T. The sample tilt was fixed 
at a value maximizing the diffracted intensity in the 'top' spots, and data were taken versus either 
increasing or decreasing temperature. The results from SANS-I were the most complete; the data 
taken at other instruments were included after scaling in the overlapping regions [H]. It is clear 
that all temperature dependencies are in agreement within statistical errors in all cases, indicating 
that the system is in the thermal equilibrium state and the vortices are not relaxing to a new state 
on temperature cycling, (b-d): Fits of the data to various gap functions versus T/T C1 with T c set 
to the value T C 2(0.1 T) = 3.25 K. The fitted parameters are shown inset in the graphs: (b) single 
full gap and single nodal gap; (c) single nodal gap including nonlocal effects [14J. These cause a 
crossover from T-linear to T 2 behavior below a temperature T* [55] . Similar effects can occur due 
to impurity scattering [31], but are not expected to be strong in our clean material. 
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